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Abstract. - The ground state tunnel splitting for the biaxial spin model in the magnetic 
field, TL — —DS X + ES Z — gfiBS z H z , has been investigated using an instanton approach. We 
find a new type of spin instanton and a new quantum interference phenomenon associated with 
it: at a certain field, Hi = 2SE 1 ' 2 (D + E) 1 ' 2 /(g/Js), the dependence of the tunneling splitting 
on the field switches from oscillations to a monotonic growth. The predictions of the theory 
can be tested in Fes molecular nanomagnets. 
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£»o I During the last years much effort, both theoretical and experimental, has been devoted to 

£N| ■ the study of macroscopic quantum tunneling (MQT) in spin systems [1, 2]. Crystals formed 

by weakly interacting identical magnetic molecules, like Mni2-acetate and an octanuclear iron 

cluster Fes, have been the objects of most intensive recent research. EPR [3] and neutron 

scattering [4] data show that an Fes cluster can be described by the biaxial spin Hamiltonian, 
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H = -DSl + ESl-9n B S-H , (1) 



i-rt ■ where S = 10, D w 0.23-0.27 K, and E s=s 0.093 K [3]. Macroscopic magnetic measurements 
[5] have revealed resonant spin tunneling in Fes, confirming the above model and values of 
the constants. At H = this model has been studied by a number of authors[6, 7, 8, 9]. The 
tunnel splitting A of the ground state in the case H ^ has been theoretically studied using 
the instanton method by Garg [10], who discovered an interesting topological effect: oscillation 
of A on the magnetic field applied along the hard magnetization axis. Such oscillations have 
been recently observed in Fes by Wernsdorfer and Sessoli [11]. The origin of this effect is in 
the quantum interference of different instanton paths, suggested earlier, within the context 
of MQT, by Loss, DiVincenzo, and Grinstein [12] and by von Delft and Henley[13]. This 
oscillatory behaviour has also been derived using standard perturbation theory by Weigert[14]. 
The model of Garg has elucidated the fact that the freezing of tunneling need not be related 
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Fig. 1 - The tunnel splitting A as a function of field h at S = f and A = O.f . Dashed line marks the 
unfreezing field hi = 0.949. 



to Kramcrs's degeneracy. To illustrate his point, Garg studied the field range < h < hi — 
(1 — A) 1 / 2 , where h — H/H c , H c = 2S(D + E)/(g/j,B) is the critical field at which the energy 
barrier disappears, and A = D/(D+E). In that field range the tunneling is quenched whenever 

[10] 



h = h 2 {S-n-l/2)/S , 



0,1, 



,25-1. 



(2) 



This was discovered before the relevant system, Fe§, was known. Meantime, for Fe 8 A w 0.71- 
0.75 and, thus, the above field range covers only the lower half of the range, < H < H c w 4.8- 
5.5 T, available for the tunneling studies. In this Letter we have considered the remaining field 
range, hi < h < 1, and found another interesting topological effect: switching from oscillations 
to the monotonic growth of the tunnel splitting on the field. We then give explanation to 
such a behavior within a continuous spin model containing the Wess-Zumino-Berry term[15]. 
Hamiltonian (1) has been recently studied by Kou et al.[16] who found exact instantons at 
h < hi but did not attempt to solve the problem in the field range hi < h < 1. We show that 
there exists a new type of spin instanton in that range, that involves motion in both imaginary 
and real time, which is ultimately responsible for the above peculiar dependence of the tunnel 
splitting on the field. Complex-time instantons have been used for some time to study barrier 
penetration effects using the path integral formalism[17, 18, 19]. To our knowledge, though, 
this is the first time that such a mixed-time trajectory is needed in the study of tunneling 
in spin systems. Garg's and our predictions can be quantitatively tested, without any fitting 
parameters, in experiments on Fcg. 

Before studying instantons we will demonstrate the reality of the effect by performing the 
numerical diagonalization of the Hamiltonian (1) for S = 10 in the field applied along the 
hard axis Z (see also rcf. [20]). Figures 1-3 show the results of numerical computations of 
the ground state tunnel splitting as a function of h for three different values of A: small A, 
intermediate A that corresponds to Fes, and large A. It can be clearly seen from the figures that 
Garg's oscillations exist below a certain field Hi. Above that field there are no oscillations; 
the tunnel splitting grows linearly with the field. The value of H 2 increases with decreasing A. 
Below, we obtain this effect via instanton method and show that the crossover field is given 
by /i 2 (A) = (1 - A) 1 / 2 . 
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Fig. 2 - The tunnel splitting A as a function of field h for Fes, S = 10 and A = 0.72. Dashed line 
marks the unfreezing field hi = 0.529. 



In the continuous approach the tunnel splitting can be computed as the functional integral 

(3) 



P[cos0(i)p[0(t)]exp( - / dtC[6(t),(f>(t)} 

over closed trajectories which describe motion of a fixed-length vector of the magnetic moment 
M in spherical coordinates (9(t), <j>(t)). Here C is the Lagrangian of the magnetic system, 



£ = — (cos(0)-l)^-W(0,« 

7 



(4) 



and 7i is the energy of the system, which includes the magnetic anisotropy and the Zecman 
term, 



w(M) 



-Mfc||sin 2 (6l)cos 2 (^) + -Mk ± cos 2 (9) 
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Fig. 3 - The tunnel splitting A as a function of field h at S — 10 and A = 0.99. Dashed line marks 
the unfreezing field h-z = 0.1. 
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1 / H 2 \ 

-MHcos(9) + -M |fc||+- — . (5) 

2 V n k± + fey y 

Here 7 = g/iB/ti is the gyromagnetic ratio, if is the external field, and fc^ and fc|| are the hard 
and easy anisotropy constants, correspondingly. The correspondence with the parameters of 
the Hamiltonian (1) is the following 

2SD , 2SE 

, k± = 

g/j-B g^B 



fc|| = , fcj_ = , H c = k± + fc|| 



* = f . A=A' («) 

fry fey + k± 

For ii < H Cl the magnetic energy Ti(0,(p) has two equivalent minima at = 0,7r with 
cos(#) = H/H c . A constant has been added to make the energy of the minima zero. The 
following observation is relevant to the calculation provided below. The energy (5), besides 
having the minima mentioned above, also has two equivalent extrema at <f> = ±7r/2 and 
cos9 — h/hi, where hi = 1 — A. At < h < hi these are the maxima of the energy. They 
would become local minima of the energy at hi < h < hi and global minima at hi < h < 1 
if one allowed for complex 0, because cos 9 = h/hi > 1 in that field range. We shall see that 
this is the case for the effective potential in the quantum problem. 

After Gaussian integration over cos(9) in eq. (3), the remaining functional integral over <p 
contains the imaginary-time effective action 



^ = ^jdTC E = -iSJdcp[A{4 > )-l] 



+SX 1 / 2 dr' 



\M(<t>)ft, + V{cj>) 



(7) 



where Le = —C(t — » —it) is the Euclidean Lagrangian, and r' = T7[fci|(fcj_ + fcn)] 1 ^ 2 is the 
dimcnsionless imaginary time. The functions A(cj>), M((f>) and V((j>) are given by 

h 

1 -Asin 2 (^) 



A (4>) = ; — T7T, - (8) 



M ^ ] = 1 An^ ' (9) 

1 — Asm ((f>) 



and 



y(0) = isin 2 (0) 



1- 



h 



j 



(10) 



1 - Asin 2 (0). 

The action (7) is equivalent to that describing the motion of a particle of mass M(<j>) in an 
inverted scalar potential —V((j>) and a "vector" potential A(<f>) — 1. The shape of the potential 
V((f>) for three different ranges of the field is shown in fig. 4. For h < hi = 1 — A, the potential 
looks like a regular barrier between = and <f> = it, with a maximum at <fi — w/2. For 
hi < h < hi = (1 — A) 1 ' 2 , the maximum at <f> = 7r/2 becomes a local minimum, though still 
higher than the global minima at <j) = 0, w. For hi < h < 1, the minimum at (p = 7r/2 becomes 
the global minimum. 

For < h < hi one can find the imaginary-time instanton trajectories for the particle 
moving from <j> = to = ±7r in the inverted potential. For these trajectories, the first 
integral in eq. (7) gives an imaginary Wess-Zumino-Berry contribution to the action, 

\SI% ZB = -isjdcp L4(0) - 1] = ±mS (l - £) • (11) 
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Fig. 4 - The potential V(<f>) at A = 0.72 (hi = 0.28 and h 2 = 0.529) for three different ranges of the 
field. The chosen fields are: (a) h = 0.2, (b) h = 0.45, and (c) h = 0.7. 



The interference of these two trajectories in the functional integral gives a factor 

h 



en I wzB -\- e k I wzB — 2 COS 



ITS 1- 



hi 



(12) 



which is responsible for the non-Kramers freezing of tunneling at fields satisfying eq. (2) [10]. 
However, at ft-2 < h < 1, one cannot find a trajectory in imaginary-time connecting and 
7r. There seems to be instead a bounce trajectory from cp = to fa and from fa to 7r, where 
sin (0i,2 ) = (1 — h 2 )/\. We then consider the instanton given by the trajectory which consists 
of three parts: 

1. motion in imaginary time from to fa, 

2. motion in real time from fa to fa, 

3. motion in imaginary time from fa to n, 

and another instanton given by the trajectory going in the opposite direction. It is clear 
from the shape of the potential that all three parts of the trajectory join smoothly, because 
far) and fa(r) coincide at the joints. Note that the real-time part of the trajectory still 
corresponds to the virtual rotation of the magnetic moment, because for that part of the 
trajectory cos(6*) > 1, as can be seen from the classical equations of motion that follow from 
the Lagrangian (4). 

Using the energy conservation 



1 



H(<f>,fa,) = -M(fafa Tl -V(fa=0, 



(13) 



one obtains 



A 2 



2V(fa 
M(fa 



= (l-h 2 )sm 2 (fa 



1 



1-h 2 



sin 2 (0) 



(14) 



This can be used to compute the second integral in eq. (7). For h^ < h < 1, the real part of 
the imaginary-time action is given by the contribution of the bounce trajectories from to fa 
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and from <f>2 to tt. Integration gives 
Re [£\ = 2SX 1 ' 2 



45 



d4>[2M(4>)V{4>)] 1/2 



where 



z = cosh 



[0,4>l]U[^2,7T] 

1 \ 1/2 ^O 2 -^^ 

— + 1 tanh - 

a 2 J I 1 + a 2 



1-A 



(15) 



h 2 - (1 - A) ' 



& 2 = 



1-A 

A 



fen 



(16) 



Equation (15) is the WKB exponent for the tunnel splitting. 
In the uniaxial limit (A— »1), one obtains 



Re I - ) = AS 



cosh 



rV(i-* a ) 1/a 



(17) 



in full agreement with ref. [21]. This formula is correct in the entire field range < h < 1. In 
the limit of small barrier[22], h = 1 — e with e^O, it gives 



Re 



8\/2 



Se 3 / 2 . 



(18) 



For h <hi-i the imaginary part of the action for the trajectory going from <fi — to <j> — ir is 
given solely by the topological term {i/K)SI^ rzB . For h 2 < h < 1, however, we have another 
contribution coming from the second integral in cq. (7) evaluated for the real-time motion 
between <pi and <j>2- Both contributions combine into 



Irn 



lCT± + C\l/2 



<i>2 



dcj)[2M{(j))V((t))} 1/2 



- -ST* 

— h WZB 



inS [ 1 - — I = 



(19) 



Remarkably, the real-time motion exactly cancels the contribution of the Wess-Zumino-Berry 
phase. Consequently, the freezing of tunneling does not occur at high fields (/12 < h < 1) and 
topological oscillations are suppressed, as confirmed by our numerical diagonalization of the 
Hamiltonian at S = 10, figs. 1-3. In each of the figs. 1-3 the field hi seem to coincide with 
the inflexion point on the envelope curve. 

Formally, one could argue that at hi < h < 1 the instantons connecting the <f> = ±7r/2 
global minima of the effective potential shown in fig. 4 should be used to compute the tun- 
neling splitting of the ground state. This would be conceptually incorrect, however, because 
tunneling in the biaxial model occurs between the energy minima <f> = 0, w, which are the 
only classical energy minima of the energy. Besides, the tunneling splitting computed via the 
(j) = ±7r/2 instantons is different from the one given above and checked through the numerical 
diagonalization of the Hamiltonian. 

One should notice that for hi < h < 1 the particle spends finite real time under the 
barrier, though the motion is still virtual because cos(#) > 1. This time is given by T = 
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27r/[7[(l — h 2 )k« (fen +fcj_)] 1 ' 2 ]. It tends to the period of small oscillations around the minimum 
at 4> = 7r/2 when h — ► hi, and diverges when the field approaches the critical value h = 1. 
Notice also that complex instantons resembling our instantons have been recently introduced 
in a rotating black hole problem [23]. 

Taking the average of the measured values of the parameters[3, 4, 5], we get A w 0.72 
for Fes, which gives hi « 0.53, and H c s=s 5.1 T. Then the Garg's oscillations must exist at 
H < Hi k, 2.7 T, while above 2.7 T the tunnel splitting must grow monotonically with the 
field. This prediction can be tested in EPR and inelastic neutron scattering experiments on 
well oriented crystals of Fe 8 . 
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